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Abstract. We present some developments in determining Ha distances 
to high-velocity clouds (HVCs) in the Galactic halo. Until recently, it 
was difficult to assess the nature and origin of HVCs because so little was 
known about them. But now several HVCs have reliable distance bounds 
derived from the stellar absorption technique, and more than a dozen 
have abundance measurements. In addition, twenty or more HVCs have 
been detected in Ha (and a few in optical forbidden lines). Over the past 
five years, we have been developing a model of the halo radiation field 
which includes contributions from the stellar disk, the stellar bulge, the 
hot corona, and the Magellanic Cloudsj^ In certain instances, the Ha flux 
from an opaque Hi cloud can be used to derive a crude distance constraint 
to the cloud. For a UV escape fraction of /esc ~ 6% perpendicular to the 
disk (/esc ~ 1 — 2% when averaged over solid angle), the HVCs appear 
to be broadly consistent with the spiral arm model. We caution that a 
larger database with full sky coverage is required before the usefulness of 
Ha distances can be fully assessed. We present a number of detailed pre- 
dictions from our distance frame to encourage independent assessments 
from future observations. If the model is valid, we find that most HVCs 
detected to date are scattered throughout the halo up to distances of 50 
kpc from the Sun. Most of this material is likely to be debris from re- 
cent galaxy interactions, or even debris dislodged from the outer Galaxy 
disk. We propose some future tests of the Ha distance model and briefly 
discuss recent Ha detections along the Magellanic Bridge and Magellanic 
Stream. 



1. Introduction 

Observations of the Galactic halo make a compelling case that the formation of 
halos continues to the present day (Wyse 1999). The halo appears to have built 



^The diskhalo ionization code, along with full documentation (Bland-Hawthorn & Maloney 
2001), is to be made available for general use. 
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up through a process of accretion and merging of low-mass structures which is 
still going on at a low level. Hierarchical cold dark matter (CDM) simulations, 
however, predict that the Galactic halo should have many more satellites than 
are actually observed (Klypin et al. 1999; Moore et al.l999). 

As much as 40% of the sky is peppered with high- velocity Hi clouds (HVCs) 
which do not conform to orderly Galactic rotation (Wakker & van Woerden 1991; 
Putman 2000). These arc interesting accretion candidates — particularly if they 
are associated with dark matter 'mini halos' — except that their distances, d, are 
unknown for all but a few sources. As a result, fundamental physical quantities 
— size (oc d) and mass (oc c/^) — arc unconstrained which has encouraged wide 
speculation as to the nature of HVCs (Wakker &: van Woerden 1997). 

Indeed, the current renaissance in HVC studies can be traced in part to 
an interesting suggestion by Blitz et al. (1999). They showed that the velocity 
centroids and groupings of positive/negative velocity clouds on the sky may be 
understood within a reference frame centered on the Local Group barycenter 
(c/. Zwaan & Briggs 2000). They interpret HVCs as gas clouds accreting onto 
the Local Group over a megaparsec sphere. Braun & Burton (2000) identify 
specific examples of compact clouds that have 'rotation curves' consistent with 
CDM mass profiles. For sources at 700 kpc, the kinematic signatures imply a 
high dark-to-visible mass ratio of 10—50. 

This paper addresses the distance problem. In earlier papers {e.g., Bland- 
Hawthorn et al. 1998; Bland-Hawthorn & Maloney 1999 - BM99), we showed 
that faint Ha measures from distant Hi clouds could, in principle, be used to 
estimate crude distances to the clouds. The Ha emission measure from any 
cloud which can be detected at 21 cm is a direct measure of the Lyc (Lyman 
continuum) radiation field, independent of distance. In order to interpret the Ha 
emission, we require a realistic model of the Galactic halo ionizing field. Here, we 
present some recent developments in deriving Ha distances which significantly 
extends our earlier work. 

2. The Escape of UV radiation from the Galaxy 

Since we have discussed the UV escape fraction (/esc) from the Galaxy at length 
in earlier articles (q.v., Bland-Hawthorn & Putman 2001), only a short discus- 
sion is given here. The likely value of /esc is still uncertain but there is increasing 
evidence that UV does radiate far from HII regions. To cite one of several recent 
demonstrations, the warm (diffuse) ionized medium in external galaxies clearly 
shows the hallmark of the HII regions even after these have been masked out of 
the Ha image {e.g., Zurita et al. 2001; Cianci 2001). 

Further evidence that UV must escape the Galaxy comes from the mea- 
sured electron density profile from halo pulsars. Manchester & Taylor (2000; see 
also Nordgren, Cordes &: Terzian 1992) have modelled this with a scale height of 
800 pc which exceeds or is comparable to the scale height of the diffuse Hi (warm 
neutral medium; Lockman 1984). Without fine-tuning, it is unlikely that the 
Reynolds layer represents a radiation-bounded medium within a co-extensive 
Hi envelope. We know that the radiation field must be soft from the weak- 
ness of HeIA5876 and non-detection of [HeII]A4686 (Reynolds & Tufte 1995). 
Furthermore, the observed weakness of [OI]A6300/Ha indicates two things: (i) 
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Figure 1. The spiral arm distribution of UV emitting sources used in 
the diskhalo ionization code where the dot shows the Sun's position. 
The model is identical to the Taylor-Cordes density distribution in the 
near field (r<10 kpc). In the far field, the spiral arms are forced to 
follow the Ortiz-Lepine model shown as continuous lines. 

the ionization fraction must be high (Reynolds 1989), (ii) all of the UV pho- 
tons produced in the disk cannot be absorbed in radiation-bounded Hii regions 
(Domgorgen k Mathis 1994). 

It is important to appreciate that the Fabry-Perot 'staring' method is so 
sensitive that we only need /esc ~ 1% for the Ha distance method to be useful. 
In fact, our preliminary analysis of Ha levels from Hi clouds with stellar distance 
brackets suggests that indeed /esc ~ 1—2% (see § 6). At distances of 300 kpc 
along the polar axis, the Galaxy field is comparable to the cosmic UV field. The 
expected Ha levels are of order 1 — 2 mR and it may be possible to reach these 
levels^ with new differential techniques {e.g., Glazebrook & Bland-Hawthorn 
2001). 

3. The importance of spiral arms 

Our early attempts to derive Ha distances made use of a smooth exponential 
distribution of ionizing sources (BM99). For HVC distances within 10 kpc of 
the plane, we need a more realistic distribution of ionizing sources than the 
exponential disk model. This would be straightforward if we knew the exact 
location of all O stars, and the precise dust distribution throughout the Galaxy. 



^1 milliRayleigh = 10^/47r phot cm"^ s"^ sr"^ = 2.4 x 10"^° erg cm"^ s"^ sr"^ at Ha. 
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vertical distance, kpc 

Figure 2. A comparison of the halo ionizing flux for different disk 
distributions (uniform emissivity, exponential and spiral) within the 
diskhalo code compared to a simple inverse square law. The vertical 
distance is measured from the center of the disk along the polar axis. 
The top three curves are in the absence of dust and converge in the 
far field limit. The lower three curves include the effects of dust where 
Ti.L= 2.8. 
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But this will not be possible until the GAIA astrometric mission flies in 2010. 
However, most studies of spirality in the Galaxy agree that the tangent points of 
the spiral arms are well defined over a wide range of methods and wavelengths, in 
particular, the distribution of pulsar dispersion measures with galactic longitude 
(Taylor & Cordes 1993). 

For this reason, the non-axisymmetric component of the diskhalo ioniza- 
tion model links its fortunes to the standard model for determining pulsar dis- 
tances. Rough distances to pulsars are deduced from the dispersion (and scatter- 
ing) measure due to warm electrons along the line of sight. Early attempts used 
a smooth distribution of electrons {e.g., Manchester & Taylor 1981) although 
Lyne, Manchester &: Taylor (1985) showed that typical distance estimates have 
random errors as large as a factor of two. After the inclusion of smooth spi- 
ral arms, Taylor & Cordes (1993) predict that most distances should be good to 
~20%. This level of accuracy is somewhat surprising when one examines face-on 
spirals in the Ultraviolet Imaging Telescope (UIT) database. But the distance 
model is largely borne out by lower limits derived from pulsar sight lines which 
show Hi in absorption. 

4. Stellar bulge contribution 

A component largely overlooked to date is the UV field arising from the Galactic 

spheroid (O'Conncll 1999). A 'UV upturn' component is now well established for 
ellipticals and SOs (Greggio & Renzini 1999; Macchetto et al. 1996; Bica et al. 
1996). Binette et al. (1994) suggest this may arise from post-AGB stars evolving 
from an old stellar population. A rough estimate is 7 x lO5O(Mbuige/lO^°M0) 
ionizing phot s~^. Mbuige is the mass of the Galactic halo bulge which is un- 
certain to a factor of two, i.e., Mbuige ~ 1 — 2 x IG^^M© (Dwek et al. 1995; 
Zhao 1996). Since much of the bulge stars lie outside most of the absorbing 
ISM, the UV flux which escapes into the halo may be non-ncgligiblc and more 
smoothly distributed compared to flux which escapes the disk. When a bulge 
component is discussed below, wc include an isothermal distribution of sources 
with a total luminosity of 7 x 10^'^ phot s^^. Wc assume that all of the UV 
escapes, although the dusty disk blocks half of the radiation. The hot corona 
component is discussed in BM99. Both sources of UV are examined briefly in 
the next section. 



5. The disk-halo ionization model 

The diskhalo model includes five basic components: the spiral arm (or expo- 
nential) disk defined by the OB star population; stellar bulge; hot coronal halo; 
LMC (and SMC); cosmic background. In addition, the resulting radiation field 
can be moderated by the presence of a plane parallel opacity law, and projected 
Ha emission measures can be derived with a variety of schemes. 

In Fig. 1, we show how the spiral arms are represented in the ionization 
model. In the near field (r < 10 kpc), the distribution is essentially identical to 
the Taylor-Cordes model. However, the spiral arms are incomplete on the far 
side of the Galaxy which required us to extend the arm coverage by 60%. This 
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was done by splicing the Taylor-Cordes model onto the basic model parameters 
of Ortiz & Lepine (1993). The spiral arms cover a total area of 100 kpc^ within 
a circle with 12 kpc radius, compared with the exponential model which extends 
out to a radius of 25 kpc (although most of the emissivity lies within a 12 kpc 
radius). 

The diskhalo manual (Bland-Hawthorn &: Maloncy 2001) describes how 
the electron density field is inverted to produce the surface density of ionizing 
photons. Fig. 2 shows how the three different disk distributions compare as we 
approach the disk along the polar axis of the Galaxy. Fig. 3 presents a cross 
section through the halo field and shows that the form of the dusty spiral is very 
different to that of the exponential disk within 10 kpc, which can also be seen 
in Fig. 2. 

An important prediction of the spiral arm and the exponential models is 
that ionizing flux seen by a cloud 'inside' the Solar Circle should be much higher 
than the flux seen by a cloud on the 'outside'. We have tried to illustrate this 
effect in Fig. 4. For this trend to be evident, a large nearby sample of clouds 
is required over a wide solid angle. One possible population which might reveal 
this effect are the intermediate velocity clouds (IVCs). In the next section, we 
show the prediction for one such cloud. Complex K. 

Fig. 4 shows how a stellar UV bulge (or, for that matter, a hot UV corona) 
can wash out some of the contrast produced by spiral arms (see also Fig. 8) . The 
presence of a UV-bright stellar bulge can remove the near-field solution which 
typically arises with the spiral arm models. This phenomenon is seen in the 
exponential disk predictions in Fig. 5. The hot galactic corona is not as effective 
as the stellar bulge at washing out the spiral structure since, even though it has 
a scale length a factor of 5 larger than the bulge, the expected flux levels are 
lower by about the same factor. 

6. Predicted Ha distances to Hi clouds 

Here we present predictions for HVCs with measured Ha fluxes. We emphasize 
that the method is only intended as a statistical constraint. Even with large 
numbers of Ha detections for each cloud complex, it may only be useful to 
within a factor of a few and even this level of confidence needs to be tested. 

We used the diskhalo code with the dusty spiral disk (rj^j^=2.8). We use 
a conservative model with a hot corona but no stellar bulge component. The 
choice of r^^ is discussed in detail elsewhere but largely arises from the WHAM 
detections of Complexes A, C, K and M (see below). Weiner et al. (2001, these 
proceedings) derive a similar value from an analysis of their own observations. 
Since galactic coordinates have not yet been published for some of the detections, 
we have averaged our calculations over the observed Hi. 

A major benefit of the spiral arm model is the order of magnitude larger 
contrast in Ha that one obtains compared to the exponential disk model. A good 
example of this comes from the Complex L sight line presented in Fig. 5. This 
cloud is found to be very bright in Ha by both the Las Campanas (Weiner et 
al. 2001) and TAURUS teams (Putman et al. 2001, in prep.). The exponential 
model fails by a large factor to explain the signal whereas this is a natural 
consequence of the spiral arm model if the complex is 1 kpc from the disk. 
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Figure 3. A cross section through the halo field produced by the 
diskhalo model for two cases: the dusty exponential disk (white con- 
tours and halftone) and the dusty spiral disk (black contours). The 
Galactic Center is at (0,0). The numbers give the ionizing flux in units 
of log(phot cm^^ s^^). 



8 



Bland-Hawthorn 



Radius about Sun = 4 kpc - 
Height above disk = 1 kpc" 




-100 100 

Galactocentric angle, & 



Figure 4. The predicted ionizing flux around a ring at 1 kpc ele- 
vation with a radius of 4 kpc centered on the Sun (6 = is towards 
the Galactic Center). The upper curves show the flux for optically 
thin cases {tli^= 0) of a uniform, exponential and spiral disk; the lower 
curves are for the dusty models (r^L= 2.8) of the same. The dusty spi- 
ral disk has been computed with and without the presence of a stellar 
UV bulge. The expected flux levels are much higher in the direction of 
the Galactic Center ('inside') compared to the Anti-Center ('outside'). 
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in which case the cloud system hes directly above a spiral arm. Independent 
support for this comes from the observation by both teams that the [Nil] /Ha 
ratios are very elevated in this system, an effect seen at low latitudes in external 
edge-on galaxies. The same temperature effect (Reynolds, Haffner &: Tufte 1999) 
is seen in the Smith Cloud which is also thought to lie within a few kpc of the 
disk (see below). 

But the contrast effect comes at a price. For most of the modelled sight lines 
in Fig. 5, there is a near-field and a far-field prediction. This is because in almost 
any direction away from the Sun, the radius vector crosses a spiral arm. Multiple 
solutions are generally avoided at latitudes higher than 20° because as we move 
away from the disk, the details of the distribution become less important. An 
example is the Smith Cloud (GCP) discussed in Bland-Hawthorn et al. (1998). 
This cloud is predicted to lie at 1 kpc or roughly 15 kpc. Two of the sight lines 
(Complex L, Co-Rotate) are multiple valued for the observed Ha. 

The WHAM detections of Complexes A, M and C are also modelled in 
Fig. 5. These clouds are particularly important as they have distance bounds 
from the stellar absorption line (SAL) method. If a star of known distance lies 
beyond the cloud, and another star of known distance falls in front, we obtain 
a distance bracket for the cloud if it is seen in absorption against the distant 
star, but not in the nearer star. Complex M (specifically. Cloud Mil) has an 
upper limit of 4 kpc on its distance (Ryans et al. 1997); no lower bound exists. 
The far-field limit is only consistent at the 1.3a level. Complex A lies between 
4 and 10 kpc from the Sun (van Woerden et al. 1999) which appears consistent 
with the far field limit. There is a suggestion that this range could be tightened 
to 8 and 10 kpc should the complex not be detected on the spectrum of PC 
0832-1-675, in which case the far field prediction is a factor of two too small. 
Complex C has a secure lower limit on its distance of 1 kpc, based upon five 
stellar probes; a much weaker limit of >6 kpc is provided by the non-detection 
of the cloud in BS 16034-0114 (Wakker 2001). Both near- and far-field solutions 
are consistent with a firm distance limit (>1 kpc); only the far-field solution is 
consistent with the weaker limit (>6 kpc). 

The WHAM team have recently published Ha observations for the IVC 
Complex K (Haffner et al. 2001). There is little absorption line data in this 
direction although the cloud does appear to lie between 0.3 and 7.7 kpc. The 
sight line to Complex K lies directly over the tangent point of a spiral arm 
which explains the broad distribution in Fig. 5. The WHAM team present three 
averaged Ha measurements for this cloud. In the absence of a characteristic Ha 
measurement for a fully mapped cloud, we take the peak Ha detection towards 
the center of a cloud as the most important for distance determination The 
brighter detections are entirely consistent with the SAL bounds, although the 
weaker detections would push the far-field limit out to 10 kpc or more. 



^This emphasizes a weakness in the Ha method — see Bland-Hawthorn & Maloney 1999b for 
a detailed discussion — in that a detailed mapping is normally required to reliably interpret 
what part of the cloud is being lit up in Ha. However, this is not always a problem. Weiner et 
al. (2001, these proceedings) and Tufte et al. (1998) find that the internal Ha dispersion in a 
significant number of HVCs can be small. 
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In summary, for HVCs with well defined distance bounds, we find (see also 
Weiner et al. 2001) that the observed Ha is roughly consistent with /esc ~ 6% 
{tll= 2.8) although the present uncertainties are about a factor of 2. Note that 
this escape value /esc is defined orthogonally to the disk plane; the solid-angle 
averaged value is /esc ~ 1—2%. 

Finally, the Las Campanas team (Weiner et al. 2001) have achieved some 
faint detections to clouds which are thought to be farther afield. We show the 
predictions for three of these (Co-Rotate, GCN and Population N). We have 
run the code for all HVCs with Ha detections to date. Population N appears to 
lie at the greatest distance (30—45 kpc) if the Ha emission arises from the disk 
radiation field. 



7. The Magellanic System 

The LMC has several highly active star forming regions, particularly regions 
of very recent star formation (Shapley III) and of ongoing star formation (30 
Doradus). The basic ionizing requirement of the LMC from combined UV, 
optical and radio studies appears to be 5 x 10^^ phot cm^^ s^^. Within a factor 
of two, this is consistent with OB star counts (Walborn 1984; Parker 1993), radio 
continuum observations (McGee, Brooks Sz Batchelor 1972; Israel & Koornneef 
1979), and vacuum ultraviolet observations (Smith et al. 1987) of the LMC. 
However, the total number of ionizing photons produced by the LMC Hii regions, 
spread over a 5 kpc region, may be as high as 1.5 — 3 x 10^^ phot s~^ (Smith 
et al. 1987). OB star counts around 30 Dor {e.g., Parker 1993) could well 
underestimate the total ionizing flux by a substantial factor. Kennicutt et al. 
(1995) suggest that fully one third of the ionizing radiation in the LMC arises 
from within 0.5° of 30 Dor. The ground-based results may suffer from crowding 
which means that the total number of stars is underestimated. 

We now test whether the Magellanic Hi Bridge can constrain the escape 
fraction of UV photons from the Magellanic Clouds. Fig. 6 shows one such 
model prediction for the Bridge. Fujimoto & Sofue (1976) give specific positions 
for the LMC and SMC in Galactic coordinates. Here, we take the total ionizing 
flux from the LMC to be 1 x 10^^ phot s^^ and the SMC to be an order of 
magnitude smaller. For both galaxies, we assume that 15% of the UV escapes. 
The expected level in the Bridge is of order 50 — 100 mR, an order of magnitude 
higher than levels produced by the Galaxy at that large polar angle. 

There are two published claims of Ha in the Magellanic Bridge. Johnson, 
Meaburn &: Osman (1982) claim to see diffuse Ha across the entire HI bridge ex- 
tending over 30° or more at an observed surface brightness of 8 R±4 R. Marcelin, 
Boulesteix & Georgelin (1985) detected the Shapley wing of the SMC in Ha at 
the level of 4 R. These Ha detections, which are partly confirmed by the UKST 
Ha survey (Zealey 2001, personal communication), are more than an order of 
magnitude stronger than expected in our simple model proposed. It seems that 
some of the bright Ha regions arise from embedded UV sources, e.g., UV-bright 
stars in the Shapley wing. However, this does not appear to be the explana- 
tion for recent detections of bright Ha along the Magellanic Stream (Weiner &i 
Williams 1996; Weiner et al. 2001, these proceedings). 
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Figure 5. Predicted emission measures along the radius vector to 8 
HVCs. The densely shaded band arises from the dusty spiral disk model 
averaged over different cloud orientations; the lighter shaded band is 
for a dusty exponential disk. The horizontal lines show the average 
and range of emission measures for the modelled position on the cloud. 
The vertical lines show where the horizontal lines intersect the model. 
For the spiral model, there is almost always a near-field and a far-field 
solution although low latitude vectors can produce multiple solutions. 
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Since the Magellanic Stream passes directly over the South Galactic Pole 
(illustrated in Fig. 7), BM99 had hoped this would constrain f^sc from the Galaxy 
but their original model underestimates the required flux by at least a factor of 
5. Weiner et al. (2001) claim detections to some points of the Stream as high as 
1 R. The Galaxy fails to produce this level of ionization by at least an order of 
magnitude. The expected Hq levels produced by the Galaxy are shown in Fig. 8 
where we also show the additional contribution from the LMC. Since the LMC 
and much of the Stream lie close to y = 0, the model predictions are given in this 
plane. Bland-Hawthorn & Putman (2001) discuss various dynamical scenarios 
which may account for the bright Stream detections. 

In summary, the Magellanic Stream poses a clear problem for the Ha dis- 
tance method since, in a few locations, it appears to be an order of magnitude 
brighter than predicted by the model. The Magellanic Bridge is also Ha bright 
in places, but this may be due to internal UV sources. However, the HVCs 
appear to be broadly consistent with the spiral arm ionization model, and we 
find that most HVCs detected to date are scattered throughout the halo up to 
distances of 50 kpc from the Sun. Our suspicion is that most of the HVCs, 
like the Magellanic Bridge and the Magellanic Stream, result from galaxy in- 
teractions, specifically dwarf galaxies disrupted by the Galaxy. An interesting 
prospect that we have been investigating with B.K. Gibson is that some of the 
HVCs have been dislodged from the outer disk due to the disruptive passage 
of a nearby dwarf. There are several observations which support this: (i) the 
metallicities are consistent with the outer disk; (ii) the outer HI disk looks highly 
disturbed (Burton & Te Lintel Hekkert 1986); (iii) in some instances the HVCs 
show continuity in velocity with the outer disk {e.g., Gibson et al. 2001); and 
(iv) the Ha distances place many of the clouds on 10 kpc scales. We have begun 
to carry out hydrodynamical simulations in order to test this idea further. 
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Figure 6. The predicted Ha emission measure along the MageUanic 
Bridge in units of log(mR). The LMC is to the left and the SMC is 
to the right. The upper curve is due to escaping radiation from the 
Magellanic Clouds. The lower curves are due the dusty exponential 
and spiral disks {tll= 2.8). 
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Figure 7. An illustration of the LMC and the dominant clouds in 
the Magellanic Stream (Mathewson & Ford 1984) projected onto the 
Galactic X-Z plane. The orbit of the Stream lies closer to the Great 
Circle whose longitude is Z = 285°. Wc ignore small projection errors 
resulting from our vantage point at the Solar Circle. The angle 5 is 
measured from the negative X axis towards the negative Z axis where 
S = -b{0° <S< 90°) and S = b+ 180° (90° <S< 180°). 
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Figure 8. The predicted Ha emission measure along the Stream 
as a function of 5 (Fig. 7) in units of log(niR). The cloud positions 
arc illustrated in Fig. 7. All curves show the ionizing influence of the 
Galaxy (tj,^ = 2.8). The short dashed curve includes the contribution 
of the LMC; the dotted curve includes the contribution of a UV-bright 
stellar bulge. The solid curve includes the effect of the LMC and a 
stellar bulge. 
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